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ABSTRACT 

Analysis of ER7S-1 multlspectral data by automatic pattern 
recognition procedures is applicable toward grappling with cur- 
rent and future resource stresses by providing a means for re- 
fining existing geologic naps. The procedures used in the current 
analysis already yield encouraging results toward the eventual 
machine recognition of extensive surface lithologic and topographic 
patterns. Automatic mapping of a series of hogbacks, strike 
valleys, and alluvial surfaces along the northwest flank of the 
San Juan Basin in Colorado can be obtained by minimal man-machine 
interaction. The determination of causes for separable spectral 
signatures is dependent upon extensive correlation of micro- and 
macro field based ground truth oi nvations and aircraft under- 
flight data with the satellite £?•*. 


1. INTRODUCTION 


The results described in this paper pertain to automatic data pro- 
cessing of ERTS-1 MSS data (Scene I.D. 1119-17204, November 19, 1972) 
for mapping physiography and geology of an area in southwestern Colorado. 
A block of approximately 200 mi 2 located peripherally to the San Juan 
Mountain Test Site, an area defined by LARS for conducting ERTS research, 
was chosen for computer analysis. This area, centered around Durango, 
Colorado was originally chosen because it appeared as a large, cloud-free 
block on otherwise cloud-covered data (Scene I.D. 1029-17195, August 
21, 1972). Subsequent analysis was performed on the cloud-free Nov- 
ember data set. Nc correlation with aircraft underflight photography 
or MSS imagery has been performed for this area. 

2. PHYSICAL SETTING OP AREA 

The area studied near Durango ranges from 6,500 feet to about 9,000 
feet MSL. Climate is relatively moderate. Kean annual precipitation 
at Durango la 19.54 Inches but diminishes with decrease in elevation. 
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Winter snow accumulation of 10 to 40 feet Is not uncommon, depending on 
altitude. Grassland or bare soil dominates below 6,000 feet, scattered 
piTfon and juniper occur from 6,000 to 7,000 feet, and pondcrosa pine and 
aspen become increasingly abundant above 7,000 feet. Orchards and row 
crops grow in the Animas River valley, and Florida Mesa southeast of 
Durango has been irrigated for pasture and crops since about 1900. 

Ceology . The San Juan Mountains are essentially an eroded, domal 
uplift covering about 10,000 mi 2 in southwestern Colorado. They are com- 
posed principally of volcanic rocks of Tertiary age, but older bodies of 
Precarabrian, Paleozoic, and Mesozoic rocks are locally exposed ir. the 
mountainous core. At many places around the margin, Mesozoic rocks occur 
as a series of tilted strata that dip towards adjacent basins, and ferm a 
scries of linear topographic highs and lows, or "hogback" ridges and 
strike valleys. The volcanic episode was succeeded during the Pleistocene 
by climatic conditions suitable for generation ar.d vigorous activity of 
numerous, large valley glaciers that have recurred during at least 3 
irregularly-spaced Intervals of relatively snort duration. Consequently, 
the geomorphic history of this rugged mountain and basin area is exceed- 
ingly but not abnormally complex. The regional erosion surface (late 
Tertiary) is in many places overlain by varying thicknesses of glacial 
till and outwash gravels. The area has been repeatedly, gently upwarped 
in recent geologic time, so that many of the glacial deposits and asso- 
ciated rocks are unstable under present slope conditions, and major 
landslides, mudflows, and rock avalanches are not uncommon anywhere in 
the San Juan region* The best general accounts of the topography and 
geology of the region are contained in Atwood and Mather (1) and Mather 
( 2 ). 


Immediately east of Durango the hogbacks and strike valleys, 
developed respectively on hard, durable sandstones and weakly resistant 
shales, are particularly prominent. The slopes developed on sandstone 
dip southeast at about 20°; thus, some strike valleys on the northwest 
side of hogback ridges are always in shadow at the sun-time of ERTS-1 
overpass. Some of the best examples of glacially alluviated benches, 
cut terraces, and floodplains of the entire region exiot on Florida 
Mesa and along the valleys of the Animas and Florida rivers. Such 
deposits are generally highly reflective to any sensor system. 

Sandstones of the hogbacks near Durango are of terrestrial origin 
and contain coal beds of fairly high quality and good calorific value 
(14,000 BTU). Mining has been sporadic, but Zapp (3) estimated measured 
reserves of 42 million short tons of bituminous coal, and total recover- 
able reserves of 1,853 million short tons (estimating 65Z recovery) 
within 3,000 feet of the present surface. 

General atraf ^raphy 2 nd topographic relations of rock units in the 
study area are shown in Table 1. 
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Tabla 1. Caucrallzcd lithologic and topographic 
tha Durango, Colorado araa. 















3. 


ANALYTIC PROCEDURE 


As seated above, the objectives of the continuing research are to 
utilize the LARSYS automatic pattern recognition procedures to obtain 
geologic maos which refine existing maps and to produce landforra simu- 
lation maps. To map geologic materials by any remote sensing technique 
it is necessary to assume: 

1. ) That subsurface materials will manifest themselves as 

spectrally separable classes at the earth's surface. 

Lithologic types, sandstone, shale, etc., are often 
obscured with a veneer of soil, vegetation, water, and man- 
made features which reflect incident sunlight to the remote 
sensor. If the subsurface lithologies are to be mipped 
based on spectral information received by the sensor, it 
must he assumed that many spectrally separable surface 
features indicate subsurface lithologic variations. 

2. ) That lithologic types are naturally segregated into a 

limited number of discrete compositional and textural 
categories which can be recognized and classified by pattern 
recognition procedures, either machine or human. This 
assumption is false, but it facilitates clustering of similar 
lithologies into discrete classes. Transitional features 
which spectrally fall between two arbitrarily chosen training 
classes will be classified as one of the two chosen classes. 
The relative spatial distribution of these "misclassif led" 
features will reflect the natural lithologic composition 
(e.g., an area of sandy shale may be classified as a random 
distribution of feature elements of the two discrete classes, 
sandstone and shale). 

Classes Considered : Based on the information provided by Zapp 

(1949) , it was decided to divide lithologies of the Durango area into 
three discrete types for analytic purposes — sandstone, shale, and 
alluvium. Sandstone is the subsurface rock type on the dip slopes of 
the hogback ridges; shale on the cret slopes and along the strike 
valleys; alluvium on the terraces and floodplains. Considered as a 
unit, these three material type: are sufficient to describe the geologic 
materials of all erosional and deposltional surfaces in the study area. 
Since ERTS-1 passes over the Durango area at about 10:30 hours, the 
west-facing cret slopes and other west-facing declivities are in shadow 
on the ERTS-1 data’ sets. Shadow areas have significantly different 
spectral signatures than sunlit areas and have been considered as a 
fourth major class. 

To obtain training field coordinate numbers for the classes of 
interest, it was necessary to generate a coordinated graphical display 
of portions of the study area. The only ERTS-1 data set available for 
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study of the San Juan test site was an August 21 frame. Four areas within 
the study region were chosen for non-supervised classif lcation (NSCLAS) 
and display. The coordinates of these areas were procurred from televi- 
sion gray scale imagery (EDIT) of the entire ERTS frame. The four chosen 
areas contained representatives of all the rock types and topographic 
forms present in the study area. The resultant NSCLAS display provided 
reliable visual separability of shadow patterns and alluvial areas. Ap- 
plication of gcomorphlc interpretation to the shadow patterns yielded 
highly probable locations for the sandstone dip slopes and shale valleys. 
Field coordinates were then obtained for the sandstone, shale, alluvial, 
and shadow areas. Training fields on each sandstone dip slope and shale 
valley were separated into dlscrece training subclasses within the major 
lithologic classes. Nine groups of training fields — four sandstone, 
three shale, one alluvium, and one shadow — were then exposed to the 
LARSYS supervised clustering algorithm (ClASS) and the entire study area 
classified. The resultant display map showed a much nigher visual sepa- 
rability of the sandstone and shale areas than the NSCLAS display from 
which the training fields were selected. Training class performance of 
the features within the training fields was 63.72. However, th results 
were still insufficient to meet our personal satisfaction criteria. With 
refinement of the coordinates of the cralnlng fields, the addition of two 
alluvial subclasses, one shale subclass, and two shadow subclasses, 
visual correlation with Zapp's geologic map improved greatly as did the 
training class performance (89.82). This classification scheme of four- 
teen subclasses in the present degree of refinement defines the state 
of the ort of the continuing analysis at the time of writing. 


4. RESULTS 

The LARSYS routines available for visual display of classified 
areas are: 1) PHOTO, a color or gray scale coded television linage," and 

2) DISPLAY, an alphanumeric coded printer image. Both modes of classi- 
fication display were choaen for visual presentation of the results 
derived trom the Durango area experiment. PHOTO Images of the study 
region a’’e shown in Figures 1 and 2. 

Figure 1 is a physiographic simulation map displaying the four 
major classes. From darkest to lightest the classes are: shadow, 
alluvium, shale, and sandstone. The three lightest classes represent 
areas 1:. lunlight and are displayed with low contrasting gray scales. 

The darkest class, shadow, is represented as a highly contrasting black 
area. This four class display produces the illusion of the third dimen- 
sion in the perception of an observer familiar with vertical or low 
oblique remote sensing images. This effect consequently gives one a 
"feel" for the physiography of the study area by forcing perception to 
selectively focus on contrasting gray scale patterns, representing 
topographic variations which are surface manifestations of subsurface 
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lithologic variotioi.9. A color 
coded PHOTO disyiay which simu- 
lates the USGS "shaded" topo- 
graphic maps is also available. 

Figure 2 i6 a machine gene- 
rated geologic map of the Durango 
area. The four major classes 
have been displayed as three 
classes which represent the geo- 
logic materials only. Based on 
geomorphic knowledge of the 
usual relationship between bed- 
rock types and their topographic 
expression, it was asc-umed that 
the areas in shadow aie underlain 
by shale. The shadow areas are 
on cret slopes that stratigraphi- 
cally fall immediately below the 
sandstone cap rocks of the 
ridges. In a region of hogbacks 
and strike valle/s deve]oped on 
sands toner and shales, it v • be 
inferred with a high deg of confidence that the surface area strati- 
graphically below the cap rock is underlain by shale. Proceeding 'n this 
assumption, we combined che shale and shadow classes as one displiyable 
unit. The resultant gray scale coded map (Fig. 2A) visually correlates 
in general form witn Zapp's map of the same area. In detail (Fig. 2B) 
the two maps become somewhat contradictory to one another. 3ased on the 
very scanty ground truth information other than Zapp's map, we have confi- 
dence that the automatically produced ma r is more reliable than the exist- 
ing geologic map. 

There are two reasons for this confidence. First, on the automati- 
cally produced map the pattern of alternating sandstones and shales 
ftriking northeaat is seen to be abruptly offaet where the Florida River 
ib superposed p 'voss the hogbacks. The cause of this offset is inferred 
to be movement of the rocks along a fault in the earth's crust. The 
rocks on the east side have been displaced to the south rnlati/e to the 
rocks on the west side of the fault. Zapp's map shows no fault in this 
area, but a zone of weaknees, such as a fault, is necessary to produce 
the conditions required frr superposition of a stream, such as the 
Florida River, across resistant hogbacks. It should be emphasized that 
Che presence of a fault is inferred from the stratigraphic patterns ae 
recognized by LARSYS automatic pattern recognition rather than fro ^ 
visual recognition of a surface lineament repreesnted on non-clas.-'fVjd 
gray scale or color imagery. 

The second reason for confidence Is the v •» high correlation be- 
tween the machine map and Atvood'a (1932) map relict erosion surfaces 



Tigure 1. Computer generated 
physiography of the Durango. 
Colorado area. Approximate scale 
1:150,000. 


478 


* % 


;; 3 5 : 5 sSfiiiai: 



> o • 
3 I o 

H 

° « 
c ai 

° £ ' 

-O 4 J O 

o cn *j 

•D -O » 

to o e *a 1 
nobc 
in <a 
u cn 


C • * 
re w #i 

JC OOfH 

< .1 a 


rH re • 
re u o 
u a o 
o 

£ •-* * 
a re •» 
re e <•< 
u ■*< •■ 
a nOrH 


V rt ■*< 

oo .e x 

g w 2 

■rl I n. 

0 >- < 

5 2 d 

1 M • o 
ft. _ « O 

"0 E o o 
< II 3 

N TJ tl 01 H 

o > *o •• 
U 3 C --I 
rH re 

9 H < « 
.■ire -h 
re i re 
U ■ u 
n n 
<u .x 

>, *j o « 
re -h re « 

6 5 2 3 




Figure 3. Physiography of the Figure 4. Geologic rap of the 

second test area, located 30 second test area. 

miles east of Durango training 

area. Approximate scale 

1:150,000. 


Atwood shows a remnant of a Tertiary erosion surface as it exists coday 
on the dip slope of the Pictured Cliff Sandstone just east of Durango. 

The boundaries of this "ancient peneplain" remnant are surprisingly 
similar to the boundaries of the Pictured Cliff dip slope as outlined 

by the computer. These two reas ns, the firsc based on geologic inference, J 

the second on correlation with previously mapped ground truth, indicate 
a high degree of accuracy and applicability of the experimental map for 
geologic purposes. 


To test the applicability of the classification scheme to regional 
areas, an approximate 200 mi 2 area centered 35 miles northeast of 
Durango was classified by using the same statistics calculated on the 
experimental training fields located two to five miles east of Durango. 
' r bis new test area was chosen because: 1) the lithologies were similar 

... those of the study area, 2) the elevations and presumably the alti- 
tude dependent vegetative cover types were different, *nd, 3) the reflec- 
tive sun angle on the hogbacks was different. The display images of thin 
area are presented in Figures 3 and 4. The physiographic and geologic 
maps of Che test region "make sense" geologically, but no ground truth 
is available to verily our results. It is apparent that extensive field 
based micro- and macro-ground truth observations and aircraft under- 
flight data of this second area are needed before any definitive statement 
concerning the reliability and accuracy of the machine produced maps 
can be demonstrated. 
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5. APPLICATION CONSIDERATIONS 


In our opinion, -:he experimental method outlined above has several 
unique application, cor improving t «*ologic mapping capabilities. The 
degree of improvement depends on: 1) the character of the terrain, 2) 

the areal dlstrlbutli » of surface features, and 3) the purpose of the 
maps required for a particular project. Additionally: 

1. ) Multiple use of the same classification technique for producing 

alternate map types (e.g., physiographic and geologic maps) can 
be obtained from the same statistics. 

2. ) Multiple map scales, ranging from 1:?'.,000 to 1:1,000,000, 

are available to users of unaltered LaRSYS displays. 

3. ) A high degree of accuracy is obtained. 

4. ) There is a regional applicability of statistics obtained from 

a "type section" of Halted areal extent. 

5. ) Compact storage in magnetic data tape libraries is possible. 

(The Durango study area statistics necessary for display of 
the maps presented in this report are stored on about 13 feet 
of magnetic tape.) 
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